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We re-examine the effect of electron-electron interactions 
on the persistent current in mesoscopic metal rings threaded 
by an Aharonov-Bohm flux. The exchange contribution to 
the current is shown to be determined by the weak localiza- 
tion correction to the polarization. We explicitly calculate 
the contribution from exchange interactions with momentum 
transfers smaller than the inverse elastic mean free path to 
the average current, and find that it has the same order of 
magnitude as the canonical current without interactions. 
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The persistent current in mesoscopic metal rings 
threaded by a static magnetic flux <f> is a striking man- 
ifestation of the quantum nature of charge transport in 
phase-coherent systems |j],|2j . At the first sight one might 
expect that this effect can be explained within a model 
of non-interacting electrons in a random potential. Sur- 
prisingly, the seminal experiment by Levy et al. [0 re- 
vealed that in the diffusive regime the average current 
is more than two orders of magnitude larger than the 
theoretically predicted current for non-interacting elec- 
trons Q . Therefore it seems inevitable to invoke electron- 
electron interactions to explain the experiment j3). This 
was first attempted by Ambegaokar and Eckern (AE) 
||, who calculated perturbatively the correction to the 
average current to first order in the screened Coulomb 
interaction. They found an average current proportional 
to Xc^evpl/ L 2 , where A c b is a dimensionless measure for 
some suitably averaged strength of the Coulomb inter- 
action at short wavelengths, vf is the Fermi velocity, £ 
is the elastic mean free path, and L is the circumfer- 
ence of the ring. The precise value of A c b is difficult to 
estimate, because it is dominated by the non-universal 
short-wavelength part of the Coulomb interaction. Us- 
ing H A c b ~ 0.06, the theoretical current is still a factor 
of 5 too small to explain the experiment. 

More recently one of us has argued that the classi- 
cal (Hartree) contribution to the persistent current is 
strongly enhanced due to the long-range nature of the 
bare Coulomb interaction [Q. Although the arguments 
put forward in Ref. have been criticized H , we believe 
that the criticism did not properly take the essentially 
non-perturbative nature of these arguments into account 
pL However, this is not the subject of this note. 

Here we would like to re-examine the exchange (Fock) 
contribution to the average current. It is well known fllCf l 



that for disordered metals singular vertex corrections in- 
volving so-called diffusons strongly enhance the effective 
exchange interaction at small wave-vectors |q| ^ 2w/£. 
This effect has been ignored by AE who focused on 
the short-wavelength part of the interaction. The pos- 
sible relevance of these diffusive vertex corrections for 
persistent currents has been pointed out some time ago 
by Beal-Monod and Montambaux However, they 

also demonstrated that there exists an overall cancella- 
tion of the leading infrared singularities. Up until now 
a quantitative calculation of the persistent current due 
to exchange interactions with small momentum transfers 
has not been performed. In this work we shall present a 
simple solution of this problem. 

It is instructive to obtain the (grand canonical) per- 
sistent current / = —cdVL((f)/d(f) from the Grassmannian 
functional integral representation of the ratio Z / Zq = 
e -(o-si )/T £ ^ ne p ar ti-ti on functions with and with- 
out interactions. Here T is the temperature. Decou- 
pling the Coulomb interaction by means of a Hubbard- 
Stratonovich field V and integrating over the Grassmann 
fields, we obtain after the usual transformations Jl2] 
Z/Z = fV{V}e- s ° !!{v \ with the effective action |3[ 
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where the factor of 2 in front of the trace is due to the 
spin degeneracy, and / q is the Fourier transform of the 



Coulomb potential 14 1. V and Go are infinite matrices in 
frequency and wave-vector space, with matrix elements 
[V]kk> = iTV k - k , and p§ 
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Here Uq is the Fourier transform of the disorder poten- 
tial (assumed to be Gaussian white noise), and £k = 
k 2 /(2m) — fj,, where m is the (effective) mass of the elec- 
trons, and /i is the chemical potential. Identifying the 
position along the circumference of the ring with the x- 
coordinate, the flux-dependence of £l((f>) is due to the 
quantization k x = ^-(n + <fi/(f)o), n = 0,±1, ... of the 
x-component of the wave-vector. 

For a calculation to first order in the RPA (random 
phase approximation) screened interaction it is sufficient 
to expand the logarithm in Eq.(nl) to second order, 
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where JV (?) = 2T£ fc [G ] fc+g , fc and 

U (q,q') = -2r^[Go] fc + g , fe ' +? '[Go]fc'fc ■ (4) 



Note that [Gojfefe' hi Eq.(||) is proportional to <5„ n ', so 
that No(q) = S n oNo(q) can be identified with the spatial 
Fourier component of the density for a given realization 
of the disorder, and Ua(q, q') = 5 n „'IIo(q, q', iuj n ) is the 
non-interacting polarization. The linear term in Eq.(|^) 
can be eliminated by means of the shift-transformation 
V q -> V q - iT-t^f^Noiq'), where f^ A is the 
inverse of the infinite matrix with elements Sqq'f^ 1 + 
110(5,5')- To calculate the persistent current it is con- 
venient to take the derivative of \a(Z / Zq) with respect 
to <j) before performing the Gaussian integration over the 
auxiliary field V. In this way we obtain / = Iq + Ih + If, 
where Iq is the non-interacting persistent current, and 
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Note that the Hartree contribution Ih involves only the 
static part = fqOq'o °^ tne screened interaction. 

Let us now focus on the average Fock current Ip in the 
diffusive regime, where the over-bar denotes the disorder 
average. To leading order in the small parameter (kpt)~ x 
(where hp is the Fermi wave- vector) we may factorize the 
average, so that we need to know the flux-derivative of 
the average polarization dUo(ci,iu))/d(/>. Note that af- 
ter averaging the polarization is diagonal in wave-vector 
space. It is easy to show p| that the result of AE || can 
be reproduced if one retains only diagram (a) in Fig^. 
For |q| <C kp and \uj\t <C 1 (where r = Ijvp is the elastic 
lifetime) this diagram yields 
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relation <jq = (AV) -1 e 2 <7o, where V is the volume of 
the system. The prime in Eq.(^j) means that the sum is 
restricted to |Q| 2tt/£. Here the ^-component of Q is 
quantized according to Q x — (n + 2(f> / <f>o) , and L v is 
the dephasing length Jul. 

Clearly, the approximation (Q) cannot be correct for 
\q\£ <C 1, because then the density vertices in FigJ^(a) 
can be dressed by singular diffuson corrections, as shown 
in Figs.|l](b)-(d). As a consequence, this regime requires 
a special treatment. This has already been noticed by 
Beal-Monod and Montambaux JpJ . Unfortunately, the 
sum of diagrams (b)-(d) in Fig.|lfcanccls to leading order 
in (kpl)^ 1 p| , so that at the first sight it seems that 
the calculation of the leading non-vanishing behavior of 
<9IIo(q, ilo) I d<j) requires rather complicated mathematical 
manipulations. However, there is a simple way to avoid 
this technical complication: As shown in Ref. G(|, cur- 
rent conservation implies that for small wave- vectors and 
frequencies the polarization is of the form 



n (q, iuj) = A~ 
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where T>{iu>) is the generalized frequency-dependent dif- 
fusion coefficient, which is related to the frequency- 
dependent average conductivity cf(iui) via |lq] 



T>{iuj) a(iuj) 



Do 0-0 

Combining Eq.([l|) with Eqs.(||) and (0), we conclude 



(11) 



(12) 



From Eqs.© and (|T|) it is now obvious that the flux- 
dependence of the average polarization can be expressed 
in terms of the weak localization correction (||) to the dy- 
namic conductivity. After a straightforward calculation 
we finally obtain 
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where A is the average level-spacing at the Fermi energy. 
gwh(i^, <t>) is the well known weak- localization correction 
to the dimensionless average conductivity a{iuj) in units 
of the corresponding Drude conductivity ctq, 



a(iuj) = <r [l + gwh(i^, 4>)} 



(8) 
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where T>q — vp£/3 is the classical diffusion coefficient. 
Note that Vq and ctq are related via the classical Einstein 



This expression, which is one of the main results of 
this work, is valid in the regime |q| ^ 2ir/£ and A <C 



£,r _1 . For 



A more sophisticated methods are 



necessary to obtain the weak localization correction to 
the polarization |l7| . For our purpose Eq.(|l3|) is suffi- 
cient, because the average of Eq.(^J) is dominated by fre- 
quencies in the range A -C \u>\ ~ E c , where E c = T>q/L 2 
is the Thouless energy. Note that for |q|£ = O(l) Eq.(|l|) 
smoothly matches the short wavelength result (f7j). In 
fact, using T>o — vf£/3, we see that for |<?wl| <§C 1 both 
expressions are identical at \q\£ — \/3- On the other 
hand, in the regime \q\£ <C 1 the long-wavelength result 
(H) is a factor of (q£)" 2 larger than Eq.^. 
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To see whether this infrared enhancement is sufficient 
to lead to a significant exchange contribution to the per- 
sistent current, we substitute Eq.([l3]) into Eq.(^J). Then 
we obtain for the long-wavelength Fock contribution to 
the average current 

/q|nn(q, iwn) 



jlong _ C^y^yV Jld<p li0 
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where the second line is valid in the experimentally rel- 
evant limit |/ q IIo(q, iui n )\ ^S> 1. Note that the current 
( |l4| ) increases with the strength of the disorder, because 
the (negative) weak-localization correction gwL m the de- 
nominator of Eq. (|l4|) becomes more and more important 
with increasing disorder, thus reducing the screening. Of 
course, our calculation is only controlled for |<7wl| ^ 1' 
so that from now on this correction will be ignored. For 
simplicity let us assume that the ring is sufficiently thin, 
such that only the motion along the circumference is dif- 
fusive. In this case the q-sum in the square brace of 
Eq.(pT|) can be replaced by a one-dimensional integral, 
which for L ^ I is independent of the ultraviolet cutoff. 
We obtain 
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For T — > and L v 3> L the summation over the Matsub- 
ara frequencies can be performed exactly, and we finally 
arrive at 



ylong _ 

1 f — 



evp 



(16) 



where L± is the transverse thickness of the ring, and 
fU) = ly sin(47r?n^/^o) 

TT ' J m 
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for < <j>/<j)o < 1/2 . (17) 



Surprisingly, this current has the same order of magni- 
tude as the non-interacting current at constant particle 
number 
regime 



In the experimentally relevant parameter 
this current is smaller than the current due 
to the short wavelength part of the Coulomb interaction 
calculated by AE || . 

In summary, we have presented a quantitative calcu- 
lation of the long-wavelength exchange contribution to 
the average persistent current in mesoscopic metal rings. 
The current has the same order of magnitude as the non- 
interacting persistent current at constant particle num- 
ber. We have also calculated the leading weak localiza- 
tion correction to the average polarization IIo(q, uS) in the 



regime |q| 27r/£ and A ;$ |w| <C t -1 , effectively taking 
singular diffuson corrections to the density vertices into 
account. 
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FIG. 1. Diagrams contributing to 9IIo(q, iui)/d<j>. Re- 
taining only (a) reproduces the result of AE. The shaded 
triangles represent the Cooperon, the solid lines are average 
non-interacting Green's functions, and the small wavy lines 
denote density vertices, (b) Vertex correction with two ad- 
ditional diffusons (represented by shaded stripes). For small 
q and ui this diagram is much larger than (a). However, to 
leading order in (/cf^) -1 (b) is cancelled by (c) and (d). Here 
the dashed lines denote the impurity average [/ q (7_q. 
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